arXiv:1506.08538vl [cs.SY] 29Jun2015 


Multi-mode Sampling Period Selection 
for Embedded Real Time Control 


Rajorshee Raha, Soumyajit Dey, Partha Pratim Chakrabarti, Pallab Dasgupta 
Department of Computer Science & Engineering 
Indian Institute of Technology Kharagpur, INDIA 
{rajorshee.raha, soumya, ppchak, pallab}@cse.iitkgp.ernet.in 


ABSTRACT 

Recent studies have shown that adaptively regulating the 
sampling rate results in significant reduction in computa¬ 
tional resources in embedded software based control. Se¬ 
lecting a uniform sampling rate for a control loop is robust, 
but overtly pessimistic for sharing processors among mul¬ 
tiple control loops. Fine grained regulation of periodicity 
achieves better resource utilization, but is hard to implement 
online in a robust way. In this paper we propose multi-mode 
sampling period selection, derived from an offline control 
theoretic analysis of the system. We report significant gains 
in computational efficiency without trading off control per¬ 
formance. 
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1. INTRODUCTION 

Embedded software-based control systems have tradition¬ 
ally been implemented by assuming fixed sampling rates and 
fixed task periods [3]. The sampling rate is derived from 
a control theoretic analysis [9] of the system in a manner 
that guarantees desired level of control performance at all 
reachable states of the system. A uniform period can be 
implemented robustly, since we can analyze the sampling 
rate of all the control tasks and then choose an appropriate 
computational infrastructure that can statically schedule a 
periodic execution of the software components of these con¬ 
trol tasks. The schedule does not change during execution 
and hence the control performance is deterministic. 

Recent studies confirm a widely accepted belief, namely 
that a uniform sampling rate is not a good choice when mul¬ 
tiple control loops share a common computing resource such 
as an Electronic control unit (ECU). These studies establish 
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that the sampling period can be regulated to achieve sig¬ 
nificant benefits in computational performance without any 
trade off in control performance. In fact, it has also been 
shown that a non-uniform scheduling strategy can balance 
the sampling rates among the control loops sharing a ECU in 
such a way that the overall control performance improves [H 

Uniform sampling rate is typically a pessimistic choice, 
since we need to choose a sampling rate that guarantees con¬ 
trol performance at all reachable states of the system. It is 
often the case that the selected rate is necessary at only spe¬ 
cific control states of the system, whereas at all other states 
a much lesser sampling rate suffices. Adaptive sampling de¬ 
rives its benefit from this fact by intelligently regulating the 
sampling period as needed to maintain the desired level of 
control performance [T]. 

Eine grained regulation of the sampling rate may theoreti¬ 
cally determine the optimal balance between computational 
efficiency and control performance but such schemes are dif¬ 
ficult to implement in practice due to non-determinism in 
timing introduced by the computational infrastructure (in¬ 
cluding message delays, execution time variations in differ¬ 
ent paths of the control software, etc). If the sampling rate 
is known a priori then it becomes possible to develop an ap¬ 
propriate schedule for all control tasks sharing a ECU with 
adequate consideration for these types of non-determinism. 

In this paper we profess the use of coarse grained reg¬ 
ulation of the sampling rate. Specihcally we propose an 
approach where for each control loop, a limited number of 
sampling rates are chosen and a control theoretic analysis 
is used to determine the switching criteria between these 
modes. Therefore, for each control loop we have a set of sam¬ 
pling states, and an automaton that captures the switching 
between these sampling states. The global sampling state 
of the system is a concatenation of the sampling states of 
the control loops in the system. Eor each global sampling 
state, the exact schedule for the control tasks is precomputed 
considering all types of non-determinism arising out of the 
execution of software tasks. The reachable global sampling 
states are chosen based on available computational band¬ 
width and the relative priorities of the control loops. 

The primary objective of this paper is to establish the 
benefit of using multiple discretely chosen sampling rates. 
In the process, we also present the following enabling con¬ 
tributions. 

1. We outline the basis for choosing the various sampling 
rates based on use case analysis. 

2. We present an analytical approach which determines 


the criteria for switching between sampling rates so 
that control performance is not hampered. 

3. We present the construction of an automaton based 
scheduler, which implements the switching between 
sampling rates of the controller. 

We present the necessary background for the work in Sec¬ 
tion O We outline the methodology of multi-mode sampling 
period selection in Section [S] Subsequently, we use the con¬ 
trol theoretic model of an Anti-lock Braking System (ABS) 
as a running example to validate our approach. 

2. BACKGROUND STUDY 

In this section, we outline the mathematical relation be¬ 
tween the sampling period of a discrete time controller and 
its control stability. Any discrete time feedback control sys¬ 
tem can be represented as shown in Fig. mm, here x{t) is 
the input to the system, e{t) is the error signal, u{t) is the 
controller output and y{t) is the plant output fed back to 
the controller using a sensor. 



Figure 1: Discrete Time Control System 

Generally, the sampling of the continuous signal is done at 
a constant rate T, which is known as the sampling period or 
interval. The sampled signal ek = e{kT), is the discretized 
signal with /c G N. In general, for any transfer function, the 
control signal output depends on n previous control signal 
output instances and m previous error signal instances [2- 
The situation may be represented as, 

Uk — —aiUk-i — CL2Uk-2 — • • • — anUk-n + ^oe/c+ 
biCk-i + b2ek-2 + ••• + & m^k — m 

The discrete signals Uk-i,Uk- 2 , •• are the delayed versions 
of Uk by sampling period T, 2T,.. respectively. In Laplace 
domain, a time delay is introduced into a signal by multi¬ 
plying its Laplace transform by the operator Let the 

Laplace domain representation of Uk and ek be U{s) and 
E(s) respectively [ 7 ]. Hence, Uk-i,Uk- 2 , -- can be repre¬ 
sented in Laplace frequency domain as 

,.. respectively and similarly for e/c_i,.. . Thus, the corre¬ 
sponding Laplace domain representation of Eq. [1] shall be, 

U(s) = -aie“^"f7(s) - a 2 e“^^"[/(s)-+ boE{s) + 

he-^^Eis) + b2e-^^^Eis) + ... 

Substituting with the discrete frequency domain op¬ 
erator z [2] and simplifying this further we get the discrete 
time transfer function of the controller C{z) as, 

+ + + ^ UJLi{z-Zj) 

z^ + aiz^-^ H- .. + ttn ° n^=i(^ - Pi) ^ 

where Zj are the zeros and pi are the poles of the transfer 
function. Behavior of any discrete time controller can be ob¬ 
served by analyzing the poles and zeros of the corresponding 
transfer function [2]. The positions of the poles and zeros 
differ for different sampling intervals (T). Correspondingly, 
the control stability of the overall system gets effected. More 
related background is provided in Appendix lAl 


3. METHODOLOGY OUTLINE 

In this section, we outline our proposed methodology of 
multi-mode sampling period selection for embedded real time 
control. The main steps of this proposed multi-mode method¬ 
ology are as follows, 

• Step I : Developing a control theoretic model of the 
corresponding system. 

• Step II : Classification of different modes based on dif¬ 
ferent control parameters and selection of best possible 
sampling rate for the corresponding operating modes. 

• Step III : Construction of a supervisory automaton for 
controlling the mode switching. 



Figure 2: Methodology Outline 


In the following sections we demonstrate the proposed 
methodology with an extensive analysis of ABS as a run¬ 
ning example. In Section [ff] we present the control model 
of the ABS. In Section (5] we divide the driving pattern of 
a vehicle into multiple modes parameterized by vehicular 
characteristics like velocity, brake pedal pressure and slip. 
For each such mode, we choose a sampling frequency which 
ensures stability guarantee of the ABS. In Section [G] we out¬ 
line our approach for guard condition selection for switching 
between different modes, and synthesize a scheduler automa¬ 
ton which may supervise the mode selection depending on 
vehicular dynamics. In Section [T] we provide experimental 
results supporting the proposed approach. 

4. CONTROL MODEL 

ABS is an automobile safety critical driver assistance sys¬ 
tem which prevents the wheels from locking and avoids un¬ 
controlled skidding. An abstract block diagram of a vehicle 
with ABS is shown in Fig. [5l For designing the vehicle we 
used a simplified quarter car model as shown in Fig. [S] 

Here, m is the mass of the quarter vehicle, 14 is lateral 
speed of the vehicle, uj is the angular speed of the wheel, 
Fat is the vehicle vertical force, Fx is the frictional force 
transmitted to the road, Mb is the braking torque, R is the 
wheel radius and Juj is the wheel inertia. Wheel slip A is 
given as, A = 1 — ojRjVx. 

The effective braking force is dependent on the frictional 
force [T] transmitted to the road which is related to Fn as, 
Fx = -^Fn, where fi is the frictional coefficient of the road 



Figure 3: 1/4 Car Forces and Torques [T] 




















Figure 4: /i — A Curve fp 
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Figure 5: ABS Overview [T] 


design a discrete PID controller for this purpose as, 


Mb = Kpe + 


/ 


edt + Ka^ 

at 


M,{z) = [K,+ ^ + 


( 7 ) 


Here Kp, Ki, Kd are the proportional, integral, and deriva¬ 
tive gain respectively of the PID controller. T is the sam¬ 
pling interval. Mb{z), E(z) are the discrete domain represen¬ 
tation of the the braking torque, Mb and the error signal, 
e = Ad — A, which is the difference between the desired slip 
(Ad) and actual slip (A). Observe that when A = Ad then 
e = 0, i.e. A = 0. Substituting this in Eq. O Mb can be 
represented as, 

Mb = [{X-l)^-mR]V^ (8) 

rC 

Hence, it is obvious from Eq. [3 & [8] that the control per¬ 
formance, i.e. stability of the ABS controller will vary for 
different values of T, 14 and A. 


surface. Thus, as evident from Eig. [H the amount of slip 
will vary depending on road conditions. 

Relationship between wheel slip and frictional coefficient 
can be approximated, fi — /(A), using a piecewise linear 
function [6] as, 

aA, A < 0.2 

-iA+| + /3, A >0.2 

where, a G [0,8] and /3 G [—0.1,0.!]. The non-linear equa¬ 
tions for designing a quarter car model can be given as, 

(5) 

“ ^) + ^\FN^^+ ^-^Mb 

Using Taylor series expansion method [iiH] for lineariz¬ 
ing a nonlinear system we obtain a linear (affine) system 
description from Eq. [5] as, 

X = Aix Ei-\- Biu^ 
y = Cix + Diu^ 

where, x^ = [I4,A], u* ^ Mb ■ Vx , y ^ [A], Ai, Bi, Cu 
Di are system input and output matrices respectively. Ei, 
I — /(x) are the affine term and function telling the valida¬ 
tion of linearizion. The formation of this state space equa¬ 
tion from the nonlinear equations are described in details 
in Appendix m The objective of ABS controller is to decel¬ 
erate the vehicle as fast as possible, while maintaining its 
steer ability by minimizing wheel slip. The main compo¬ 
nents of ABS are the ABS-ECU, hydraulic modulator, and 
wheel speed sensor. The ECU constantly monitors the wheel 
rotational speed through the wheel speed sensors, and also 
measures the actual slip. The controller’s task is to main¬ 
tain the braking torque within a certain range. The braking 
force is applied to the wheels by the hydraulic modulator. 
It rapidly pulses the brakes to prevent wheel lock up, even 
during panic braking in extreme conditions and promises 
shortest possible distance under most conditions. We can 


5. MODE AND PERIOD SELECTION 

The candidate sampling modes and periods for a controller 
are determined by partitioning its input space based on the 
use-case scenarios and stability of the control law in those 
scenarios. Eor example, in ABS, the adequacy of a sampling 
rate in a given scenario depends on the urgency of braking 
(which is a function of vehicle speed and pedal pressure) 
and the slip ratio (which is a function of the vehicle speed 
and friction on the road). In general, we select the differ¬ 
ent possible sampling frequencies of the controller using the 
following approach. 

1. Identify vehicular parameters which impact the con¬ 
troller output (i.e. braking torque in this case). 

2. Identify multiple possible driving scenarios and corre¬ 
sponding ranges of vehicular parameters and also the 
probable driver response. 

3. Perform stability analysis followed by identification of 
maximal acceptable sampling period in each driving 
scenario. 

Steps 2 and 3 may have to be iterated to arrive at a gain¬ 
ful combination of sampling modes. A sampling mode is ef¬ 
fective towards gaining computational efficiency only if the 
controller stays in that mode for a non-trivial period of time. 
Therefore it is necessary to relate the sampling modes with 
different use-case scenarios. In ABS, we estimate different 
possible traffic scenarios (city traffic, suburban or medium 
traffic & highway traffic) and the variation in traffic den¬ 
sity, traffic regulations and corresponding average cruising 
speed and driver reaction to arrive at the sampling modes. 
We categorize the brake pedal pressure range as low, mild, 
medium and high, and also consider various speed ranges 
and slip ratios. Eor each of the scenarios, we determine the 
sampling rate at which the controller is stable. Through this 
study we selected three sampling modes, as outlined below: 
• NO Mode: This mode requires the least sampling rate 
among the three chosen mode, and targets scenarios where 
the vehicle is cruising at low to medium speeds (such as 
in city traffic). Considering average cruising speed and the 
driver reaction, we arrive at the operating sampling rate 















































Ts = 0.2ms in which the ABS achieves satisfactory control 
performance as given by Fig. [6l For a given velocity (X axis) 
and slip (Y axis), we carry out the standard unit circle anal¬ 
ysis and plot the maximum magnitude among the different 
pole positions (Z axis) of the transfer function correspond¬ 
ing to Ts = 0.2ms. The stability variation for different slip 
values is due to the piecewise linear function given in Eq. 
[4] It may be observed from Fig. [6l all the poles are of 
magnitude <= 1 for velocity range [0 .. .Sb]km/h and slip 
range [0... 0.65] thus ensuring stable vehicular dynamics. 
Thus correspondingly brake pedal pressure variation range 
is selected as [low, mild, medium] estimating the probabil¬ 
ity of the above mentioned slip range. We carry out similar 
analysis for the other cruising modes and derive satisfactory 
sampling intervals. 



Figure 6: NO Mode Stability Guarantee 

• N1 Mode: This mode uses higher sampling rate than 
NO, and targets suburban traffic scenarios. The maximum 
sampling interval with stability guarantee is found to be 
Ts = 0.15ms as shown in Fig. [3 



Figure 7: N1 Mode Stability Guarantee 

• E Mode: This mode uses a sampling rate that is adequate 
in all scenarios. Existing approaches for choosing a uniform 
sampling mode will choose this sampling rate. We choose 
the sampling frequency for which the vehicle remains stable 
considering all velocity and brake pedal pressure variations 
as shown in Eig. [HI The corresponding sampling period for 
our model is found to be Ts = 0.1ms. We designate this as 
emergency sampling mode, in case of any driving irregularity 
the controller switches to this mode, thus ensuring vehicular 
stability. 

Z-domain unit circle stability analysis is used to show the 
vehicular parameters and stability relationship graphically 
for extensive parameter ranges. Similar observations can be 
obtained mathematically using other nonlinear system sta¬ 
bility [2] criteria like, Lyapunov, Nyquist, Routh-Hurwitz or 
Bode plot analysis. The effective braking pressure range is 
varied in each of these modes to achieve satisfactory perfor¬ 
mance, depending upon driving scenarios. The mode switch¬ 
ing and detailed switching criterion are discussed in the next 
section. 



Figure 8: E Mode Stability Guarantee 


6. SUPERVISORY AUTOMATA 

Our analysis of the different possible driving scenarios and 
the choice of sampling periods entails the creation of a sched¬ 
uler which may dynamically switch the controller among 
different sampling modes. Such a supervisory automaton is 
shown in Fig. (9] The criteria for switching between the three 



Figure 9: Scheduler Automata 

modes of the automaton is chosen based on our observations 
about vehicular parameters vis-a-vis stability. Let “u” and 
be the velocity and brake pedal pressure respectively. 
The brake pedal pressure is divided in to the ranges low, 
mild, medium and high. High brake pedal pressure repre¬ 
sents the probability of larger value of slip. Medium, mild 
and low brake pedal pressure signifies lesser value of slip. 
The guard conditions for switching between states are given 
as. 
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The guard conditions for switching between WO’ to Wl’ 
and vice versa are selected to be Vni’ and V^o’ respectively. 
Similarly, for Wl’ to W’ and vice versa, switching condi¬ 
tions are Ve’ and respectively. The guard conditions 
are selected ensuring some amount of hysteresis while mode 
switching. For example, the automaton switches from mode 
WO’ to mode Wl’ in case bpp = high and v G [0...85]. 
However, the automaton switches from mode Wl’ to mode 
WO’ when v G [0... 80] and bpp is not high. Similarly, 
the automaton switches from mode ‘Wl” to E when bpp is 
medium or high and u G [> 140], however, the automaton 
switches from mode ‘E’ to mode ‘Al’ when u G [> 135] and 
bpp is low or mild. 

Whenever the automaton makes a transition from a mode 
with lower sampling period to a mode with higher sampling 
period (e.g. E to Nl), it ensures that the guard conditions 
are valid for a certain prefixed number of clock cycles. In 
that way, unwanted glitches due to faulty sensor readings 
are expected to be filtered out. 



























































































































































































































































































The main objective of synthesizing the scheduler automa¬ 
ton was to reduce ECU bandwidth requirement. Scheduling 
in E mode signifies a sampling periodicity of 0.1ms, while 
the sampling periodicity of NO and N1 modes are 0.2ms and 
0.15ms respectively. If we notice the mode switching sce¬ 
narios, we observe that when the car is cruising at a certain 
speed and no brake pedal pressure is applied, the controller 
is scheduled using infrequent sampling periods (NO or Nl). 
Further, in scenarios when the car is cruising at a certain 
speed and brake pedal pressure is applied, the mode switch¬ 
ing will be supervised by the respective scheduler automa¬ 
ton ensuring that it will switch to a more frequent sampling 
mode. Thereby, in a general cruising scenario, a multi-mode 
controller ensures nearly 30% — 50% ECU bandwidth saving 
as shown in our simulation results. 


7. RESULTS 

The initial part of this section is devoted towards estab¬ 
lishing the motivation for multi-mode sampling through ex¬ 
perimental results. The latter part of the section demon¬ 
strates the gain in computational bandwidth and the ben¬ 
efit of effective sharing of computational resources between 
multiple controllers. 

We analyze the performance of the ABS for different sam¬ 
pling intervals. Observing the // vs A curve in Fig. [4| care¬ 
fully, we notice that the peak point on most of the road sce¬ 
narios belong to the range [0, 0.2]. Hence, for effective brak¬ 
ing with maximum possible road friction, we set Xd = 0.2 as 
the desired slip. In our experimental setup, when brake force 
is applied with current velocity V = lOOkm/h, it is expected 
to gradually decrease until V = 0km/h and throughout this 
deceleration phase the slip value should be as close as possi¬ 
ble to the desired slip (A^ =0.2) thus ensuring smooth brak¬ 
ing. Finally the slip value should be I (normalized slip) 
when the vehicle comes to rest. 





Figure 10: Slip Variation: the Left &; Right column 
figures correspond to sampling time Tsi = Is sam¬ 
pling time Tsi = O.OIs respectively. 


the stability of the system. We empirically observe the vari¬ 
ation of stability with sampling rate. For stability analysis, 
we used unit circle as well as bode plot analysis [2]. We 
take a moderate velocity, say V = 60km/h and calculate 
the stability of the system for different sampling intervals. 
We show one such example Bode plot for stability analysis 
in Fig. nil It may be observed that with the sampling in¬ 
terval Ts = O.OIs, the system is unstable. However, if we 
drastically reduce the sampling interval to O.Ims, the sys¬ 
tem becomes stable. In Appendix 1C.II and 1C.21 more such 
results are given. 


Bode Diagram 
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Figure 11: Stability Variation:Ts = O.OIs, Unstable & 
Ts = O.Ims, Stable 

We have employed our methodology of multi-mode sam¬ 
pling period selection for the ABS running example. We 
consider a braking scenario with the initial and final speed 
being 200km/h and 0km/h respectively and compare the 
minimum stopping distance achieved by our ‘Multi-mode’ 
ABS Controller with the existing Matlab model of ABS with 
fixed periodicity. The simulation results are provided in Ta¬ 
ble [1] considering different possible road surfaces. We ob- 


Table 1: Stopping Distance in Kilometer 


ABS 

Controller 

Road Surface 

Dry 

Asphalt 

Gravel 

Loose 

Gravel 

Wet 

Existing 

model 

3.073 

3.424 

3.771 

4.269 

Multi-mode 

3.080 

3.433 

3.795 

4.289 


serve that the stopping distance is nearly same for both the 
controllers. For this ‘panic braking’ scenario, the estimated 
percentage of time spent in each of E, Nl and iVO modes is 
highlighted in Fig. [121 If is evident from Fig. [12] that we 
can save significant amount of ECU bandwidth (30% - 50%) 
using a multi-mode controller as compared to the controller 
with fixed periodicity. 


We observe that the slip varies significantly for two differ¬ 
ent sampling rates as given by Fig. HOI Observe from the left 
column in Fig. [10] that with a choice of moderate sampling 
rate, the slip varies drastically thus leading to undesired 
perturbations in vehicular speed while braking. However, 
as shown in the right column of the same figure, with the 
higher sampling rate, the slip exhibits a well damped trajec¬ 
tory around the desired value {Xd =0.2) leading to smoother 
vehicular deceleration. 

Our notion of control performance is based on ensuring 
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Further, we investigate the utility of our multi-mode ABS 
controller in a general cruising scenario where the car is be¬ 
ing driven in a speed range of O/cm/Zi to 200A;m//i, in various 
traffic densities. The ECU bandwidth requirement for this 
scenario is shown in Fig. 1131 It is evident that the multi- 
mode controller requires much lesser bandwidth compared 
to the existing controller since the supervisory automaton 
schedules the controller in the infrequent sampling modes 
for most of the time as described in Section [6l In that way, 
we can guarantee a significant amount of bandwidth saving 
which may be utilized for scheduling other tasks. 



Figure 13: ECU, General Cruising Scenario 

We further demonstrate the promise of multi-mode sam¬ 
pling for a multiple ECU scenario taking an Adaptive Cruise 
Control (ACC) system as a running example. ACC is an au¬ 
tomobile safety critical driver assistance system which au¬ 
tomatically adjusts the vehicle speed in order to maintain a 
safe distance from vehicles ahead. In case of an ACC sys¬ 
tem, the driver sets a safe cruising speed and a desired safe 
distance (from preceding vehicle) as controller inputs. The 
other inputs of an ACC controller which comes from the 
radar sensor are preceding vehicle speed and approximate 
distance from preceding vehicle as shown in Fig. [Tffi 
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Figure 14: ACC System Overview [T] 


ACC is a drive assist system for highway cruising which 
monitors the inputs and decides cruising speed or distance to 
lead vehicle. When the applied brake pedal pressure is high, 
the ACC is overridden by the braking controller (ABS). The 
operating modes of ACC controller can be classified into ‘ac¬ 
tive’, ‘suspended’ and ‘idle’ while our multi-mode ABS oper¬ 
ates in three different modes as discussed previously. Gener¬ 
ally, the ACC system and the braking controller are mapped 
to separate ECUs. We can achieve significant reduction in 
bandwidth requirement by sharing an ECU between these 
features. 

When the ACC is suspended, because the applied brake 
pedal pressure is high, we schedule the ABS controller with 
frequent sampling rate and the ACC controller with rela¬ 
tively infrequent sampling rates. On the other hand, when 
the ACC is active, i.e. applied brake pedal pressure is low 
or null, then the ACC controller is scheduled with frequent 
sampling rate and the ABS controller is scheduled with rel¬ 
atively infrequent sampling rate. In Fig. [T5l we show the 
bandwidth requirement of both ACC and ABS controllers 
when scheduled in a single ECU while providing satisfac¬ 
tory control performance. 
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Figure 15: ECU Sharing: ACC-ABS 


8. CONCLUSIONS 

The present work provides a methodology for adaptively 
regulating the sampling rate of embedded software based 
controllers leading to significant reduction of computational 
resource requirement. Applying the methodology on single 
controller based systems like ABS has shown that 30% - 50% 
reduction in ECU bandwidth requirement is possible. Fur¬ 
ther, it was also shown that the method smoothly scales up 
for multiple controller based systems. Our future research 
shall focus on giving a sound formal underpinning to the 
method of creating multiple sampling modes for software 
based controllers and creating a tool flow which mechanizes 
the synthesis of such multi-mode controllers. 
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APPENDIX 

A. BASIC CONTROL SYSTEM 

Any continuous feedback control system [7] can be repre¬ 
sented as shown in Fig. 1161 where x{t) is the input to the 
system, e{t) is the error signal, u{t) is the controller output 
and y{t) is the plant output fed back to the controller us¬ 
ing a sensor(H). Correspondingly, a Discrete time Feedback 



Figure 16: Continuous Time Control System 

Control System can be represented as shown in Fig. El 
where the dotted box highlights the discretized controller 
[2] . Generally, the sampling of the continuous signal is done 



Figure 17: Discrete Time Control System 


at a constant rate T, which is known as the sampling pe¬ 
riod or interval. The sampled signal yk = y{kT) is the dis¬ 
cretized signal with k E N. For understanding the design 
of a discretized control system, let us first consider a simple 
continuous domain transfer function for the controller C in 
Fig. 1161 given as follows. 


C(s) 


_ £(£) _ 
- E(s) - 


K(s-\-a) 


The corresponding time domain representation shall be, 
^+bu = K{^ + ae) 


The Euler’s approximation of the first order derivative is 
represented as: 


dx ^ Xk-\-i Xk /Q\ 

~dt ^ At ^ ^ 

Applying Euler’s approximation of first order derivative [2] 
on the continuous differential equation, we get the following 
discrete difference equation. 


+ buu = + aeu) 

Generally, At, K, a and b are fixed. The digital controller 
updates the control signal every cycle as per the following 
equation, 

Uk+i — —aiUk + h^Ck+i + hiCk ( 10 ) 

where, ho — K^bi— A"(aAt—1), ai = (5At—1). In general, 
for any transfer function, the control signal output depends 
on n previous control signal output instances and m previous 
error signal instances which can be represented as [2], 


Uk = —CilUk-l — Ci2Uk-2 — • • • — anUk-n + ^Oe/c + 
blCk-l + h2Ck-2 + • • • + bm^k — m 


( 11 ) 


The discrete signals Uk-i,Uk- 2 , • • • are the delayed versions 
of Uk by sampling period T, 2T,... respectively. 


A.l Stability: Discrete Time Control System 

The stability property of this system can be defined from 
the impulse response of a system as 

• Asymptotic stable system: The steady state impulse 
response is zero. 

lim y(k) = 0 

k^oo 

• Marginally stable system: The steady state impulse re¬ 
sponse is different from zero, but limited. 

lim 0 < y{k) < oo 


• Unstable system: The steady state impulse response is 

unlimited. .. 

iim y[k) = oo 


where y{k) is the impulse response of the corresponding sys¬ 
tem. The impulse response for different stability property is 
illustrated in Eig. El Let us assume a control system with 


stable System 




Figure 18: Impulse Response and Stability 

input u and output y. The transfer function of any discrete 
time control system can be represented as 

C{z) = 44 = 

where p is the pole which is in general a complex number 
and can be written in polar from as p = where m is 

the magnitude and 0 is the phase. The impulse response of 
the system can be given as 

yik) = Z-^{^} = b\m\^e^^o 

Thus, it is the magnitude m which determines if the steady 
state impulse response converges towards zero or not. The 
relationship between stability and pole placement can be 
stated as follows. 

• Asymptotic stable system: All poles lie inside (none is 
on) the unit circle, or what is the same: all poles have 
magnitude less than 1. 

• Marginally stable system: One or more poles but no 
multiple poles are on the unit circle. 

• Unstable system: At least one pole is outside the unit 
circle. 

The situation is graphically shown in Fig. 1191 
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Figure 19: Unit Circle: Stability areas in Complex 
Plane 

B. QUARTER VEHICLE MODELING 

A quarter car model is shown in figure [20l Here, m is the 



Figure 20: 1/4 car forces and torques 


mass of the quarter vehicle, Vx is lateral speed of the vehicle, 
cj is the angular speed of the wheel, Fjv is the vehicle vertical 
force, Fx is tire frictional force, Mb is the braking torque, R 
is the wheel radius and J^j is the wheel inertia. Wheel slip 
A is represented as: 


A = 


ujR 

Vx 


The relationship between Fn and iA[I] is given as: 


Fx = —fi{X)FN 


where /i(A) is the frictional coefficient of the road surface. 
The non-linear equations for designing a quarter car model 
can be given as [8]: 


V = -^Fn^i{X) 


A = 


_ Vx(l-X)-coR 

^ ^ - A) + ^]FNfJ^{X) + ■ 


Mb 


For linearizing the system approximation using the Taylor 
series expansion can be expressed as [6]: /(A, Vx) ~ /(A', V/) 
+ (A - A') + “ V)- From this from 

this we can derive a linear (affine) system description as: 


X — Aix + El + Biu 
y = Cix + Diu* 

I = f{x) 


wheel slip and frictional coefficient can be approximated by 
using piecewise linear function as. 


m(A) = 


aX, 

— |A + I + /3, 


A < 0.2 
A > 0.2 


( 12 ) 


where, a G [4.8,5.1,5.46,6.4] and jS G [—0.1,0.!]. 

-Ftv 


A, = 


0 

- V'J^ 


El = 


0 

.FNR^y 


and for A > 0.2 


Ai = 


(-^ + 


3 \ F]\jR 

4 A VF Joo 


± 0.1 




Fn 

Am 

FnR^ 

4V' Ju, 


El = 


(-!±o.i)± 

±0.2 


vsJu, 


(13) 


(14) 


(15) 


(16) 


The stability analysis of this system can be done using Lya¬ 
punov, Bode, Nyquist, Unit Circle or Hurwitz stability cri¬ 
teria [8]. 


C. STABILITY V/S SAMPLING PERIOD 

In this section we provide few examples relating to the 
variation in stability with respect to sampling periods. 

C.l Simple Examples: 

Let the transfer function of any arbitrary system be, 

Uis) = 

X A Yns^ -\-bs-\-u 

Case 1: Let m=2, b=-0.5, u=l. For Ts—2s the system 
unstable but for Ts=ls the system is stable. The values of 
m, b, u are same for both of the sampling intervals. The 
corresponding stability response is shown in Fig. \2T\ 


Bode Diagram 



Frequency (rad/s) 


Figure 21: Ts = 2s, Unstable; Tg = Is, Stable 


where, A/, Bi, Ci^ Di are system input and output ma¬ 
trices respectively. Ei are the affine terms and f{x) is the 
function telling the validity of linearizion and x^ = [14, A] , 
u* = Mb‘Vx — [A]. From Fig. |4]the relationship between 


Case 2: Let m=5, b=l, u=10. In this case the system 
stable for both the sampling intervals Ts—2s and Ts = ls. The 
values of m, b, u are same for both of the sampling intervals. 
The corresponding stability response is shown in Fig. [22l 












































Bode Diagram 



C.2 ABS example 

We provide more examples of stabiliy variation with sam¬ 
pling interval for the ABS controller model. 

Case 1: For V = 10 Km/h, A = 0.1, the stability response 
is shown in Fig. 


Bode Diagram 



Figure 23: Ts = 0.6 ms, Stable 
Case 2: For V = 100 Km/h, A = 0.6, the stability re¬ 
sponse is shown in Fig. 1241 


Bode Diagram 



Frequency (rad/s) 

Figure 24: Tg = 0.6 ms, Unstable; Tg =0.1 ms. Stable 


Case 3: For V = 40 Km/h, A = 0.2, the stability response 
is shown in Fig. 1251 


Bode Diagram 



Figure 25: Tg = 0.6 ms. Unstable; Tg = 0.3 ms. Stable 

Case 4: For V = 15 Km/h, A = 0.3, the stability response 
is shown in Fig. 1261 


Bode Diagram 



Figure 26: Tg = 0.6 ms, and Tg — 0.1 ms. Both Stable 

From these examples we can observe that sampling period 
has a major role to play in deciding the stability of a software 
based controller. 




















































































